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1. Introduction 

We study scalar QED with two flavors of complex, massive fields (j) and x on the lattice, where 
the masses of the two fields are chosen differently. For each flavor a finite chemical potential cou¬ 
pling to the conserved f/(l) charge is introduced to investigate condensation. In the conventional 
representation the theory has a sign problem at finite density, which we can overcome by rewriting 
the partition sum in terms of dual variables. For the theory we are studying here, a detailed treat¬ 
ment of the dual reformulation and the description of suitable simulation algorithms can be found 
in [1]. Similar techniques can be used to solve the sign problem of other theories, see, e.g., [2]. 

In this study, we are mainly interested in the connection of the condensation thresholds to the 
masses of the lowest-lying excitations of the system. This is motivated by the question whether one 
can see sequential condensation when the chemical potential is coupled to charges carried by more 
than one field, or having different quantum numbers. A particularly intriguing instance of such a 
situation is that of QCD with isospin chemical potential: When the chemical potential is increased, 
can one observe sequential condensation of the charged pion, then the charged rho, et cetera? We 
study this situation in a simple two-flavor toy model. For a related study in a theory without sign 
problem see [3]. 

Against our expectations we only find a single condensation threshold, where the original 
motivation for this study, as already pointed out, was to possibly find a sequence of thresholds. 
Nonetheless, it turns out that understanding the non-trivial relation of the mass spectrum to the 
determined critical chemical potential is interesting in its own right and we will present an adequate 
explanation for the obtained results. 

2. Conventional lattice action 

We write the action of two-flavor scalar QED with non-degenerate quark masses as a sum, 

S[u,^,z] = %[C] + S^[^,U] + S^ixM . (2.1) 

Here Su is the pure gauge action, while Sip and S), denote the actions of the matter fields (p and X 
repectively. The gauge part Sp of the action is given by the usual Wilson plaquette action, 

%[t/] = > (2-2) 

JC (7<T 

where j3 is the inverse gauge coupling, and the sum runs over the real part of all unique plaquettes 
Re Ux.aUx+d,-z^x+x ^ith Ux^a being a link variable at lattice site x in direction a. The gauge 
fields Ux.a are elements of [/(I). Here x specifies lattice sites on a 4-dimensional space-time lattice 
with extent Ns in the three spatial directions, and Nt in the Euclidean time direction. 

The matter part of the action involving the complex valued fields 0 G C is given by 

+ + . (2.3) 

X ^ v=l ^ 

We use the abbreviation = 8 -|- m^, with being the bare mass parameter of the field (p. The 
second parameter Xp is the coupling for a quartic self-interaction. Note that a chemical potential 
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/X 0 is introduced in the usual way by adding a term proportional to Q to the Hamiltonian of the 
system, where Q is the conserved charge, connected to the 17(1) symmetry of the field (/). As usual, 
the chemical potential /i 0 induces an imbalance between hopping terms in positive and negative 
imaginary time direction, as can be seen from equation (2.3) and thus the action S'o is complex 
for 7 ^ 0. The sum (2-3) runs over all lattice sites x and for the hopping terms there is an 
additional sum over all directions v = 1 , 

The action S), for the second flavor, is almost identical to S,j,, 




+ ^+Mz5v,4 


JC-V,V A.X-V 


Xx 


v=l 


(2.4) 


where = 8 + m^, and is the quaitic coupling parameter for the field X- Also here we in¬ 
troduce a chemical potential coupling to the temporal hopping terms. Note that compared to 
the hopping terms in (2.3), for the field X the link variables Ux,v enter with an additional complex 
conjugation, implying that the fields (p and X carry opposite charges. Throughout this study we use 
an iso-spin like chemical potential by setting = ix. 

The fact that for finite chemical potential and/or fx^,, this theory has a complex sign problem, 
implies that in the conventional representation the Boltzmann factor can no longer be used as weight 
in a Monte Carlo simulation. This is why we here use a reformulation of the theory in terms of 
dual variables [ 1 ], which is briefly discussed in the next section. 


3. Dual representation of the partition sum 

To avoid the sign problem of the conventional formulation as introduced in the previous sec¬ 
tion, we reformulate the theory in terms of dual variables. In the dual formulation the sign problem 
is gone and the real and positive weight factors occumng in the dual partition sum can be used as 
weights in a Monte Carlo simulation. We will here just give the results of the dual reformulation 
and for a detailed treatment refer the reader to [1]. 

In a very compact notation, the partition sum in the dual formulation can be written as 

( oo\/ oo oo oo oo\ 

n E n E E E E 

x.v<T:px„=-°°/ ;,-_v=-<»/(,v=0A,v=-°°/fv=0/ 

where the dual, integer degrees of freedom are the plaquette variables Px^vr, corresponding to the 
original gauge degrees of freedom U, the link variables I and I' represent the original d.o.f. of the 
field (p, and the link variables j and /, correspond to %■ 

From integrating out the original degrees of freedom we obtain constraints on admissible con¬ 
figurations of the system, which in (3.1) are denoted by the symbols '^g[p,l,j], ‘^s[l] and 'Tas[j]. The 
constraint "l^glp, I,j] requires closed surfaces of plaquette flux p, which have to be saturated by link 
fluxes 1 and j along open boundaries, while the constraints ^i[/] and ^s[j] demand a conservation of 
I and j flux at each lattice site. The primed link variables I' and / are not subject to any constraint. 

Admissible configurations come with a total weight >^[ 7 ,/], where, as indi¬ 

cated by the notation, the total weight factorizes into a weight Wp [p] coming from the plaquettes, 
another weight coming from the I and I' links, and a third weight >^[ 7 ,/] which depends 
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on the configuration of j and / links. In contrast to the conventional form of this theory, in the dual 
representation the weights are real and positive, even at finite chemical potentials /r^, Using the 
dual representation allows one to perform a Monte Carlo simulation of the system at finite density. 


4. Observables 


In the dual like in the conventional representation, thermodynamic observables can be obtained 
as derivatives of the logarithm of the partition sum with respect to the various couplings. The 
plaquette expectation value {U) and the corresponding susceptibility Xu^ are given by 


{U) = 


6N^Nt dp 


InZ 


Xu = 




GNfNt dp- 


InZ 


i.e., they are obtained as derivatives with respect to the inverse coupling p. The field expecfafion 
values (0*^) and (z*Z)> togefher wifh fhe corresponding suscepfibilifies X(i>*ii) and Xx*X’ 




d^ 


InZ, {x*x) = ^i^^riZ, Xrx = ^ 


d^ 


V dMl 


V d{MlY 


InZ, 


U d[Mlf 

and 

particle number densities and n^, fogefher wifh fhe associafed suscepfibilifies Xn^ and are 


and can be obfained as derivafives wifh respecf fo fhe squared mass paramefers and The 


\ d I d^ \ d \ d^ 

= C 3TU = 77 tut = 77 TTU = 77 


V dj^^ 


vd^i 


Vd^, 


V d^^ 


To defermine fhe /i = 0 mass specfrum of fhe fheory we use fhe variafional approach, where 
we diagonalize fhe connected correlation mafrix Cij{Ant) by solving fhe eigenvalue equafion 




(4.1) 


wifh fhe correlafion mafrix given elemenf-wise by 

Qj{An,) = {di{An,)d*m - (0|0,-|0)(0|0y|0) = £(0|O,-|Z)(Z|(9)|0)e-^'^' , (4.2) 

i 


and fhe zero-momenfum inferpolafors (5, defined as 

0\{nt) = ^Y^<^*{n,nt)<^{n,nt) , d2{nt) = ^Y^<^{n,nt)x{n,nt) , (4.3) 

n n 

Osirit) = ^J^(t>*{n,n,)x*{n,n,) , 04{nt) = ^'^X*{n,nt)x{n,nt) . (4.4) 

H n 

Here fhe sums run over all lattice sites H and we used fhe unconvenfional combinafions ^x 
^*X* iri the inferpolafors because fhe action given in Section 2 implies fhaf fhe fields (p and X are 
charge conjugafe fo each ofher. If can be shown fhaf fhe eigenvalues are related fo fhe masses 
El of fhe physical slates by 

A^^^(An,) oc exp (—AniA"/)(l + i^(exp(—A/ijAE;)) , (4.5) 
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Figure 1: We show results for the relevant eigenvalue and the corresponding eigenvector for the 
parameters j3 = 0.7, A = 1, = 5.3, = 5.7 at xNt = 16^ x 20. In the l.h.s. plot we performed a fit 

to the data points at An, = 1,2,3. 


and further that the corresponding eigenvectors encode the contributions from the interpolators 
in the correlation matrix Cij to the physical states, 

. (4.6) 

i 

Note that the used interpolators correspond to two-particle bound states, since we probe the system 
in the confined phase and expect the lowest excitations to be of mesonic type. 


5. Results 

On the l.h.s. of Figure 1 we show results for the eigenvalue A^^^ of the connected correlation 
matrix {Cij} given in (4.2), while on the r.h.s. we show the corresponding eigenvector as 
function of An,. The used parameters were j8 = 0.7, A = 1, = 5.3 and = 5.7 at a lattice size 

0fiV3 

xNt = 16^ X 20. We used 10^ sweeps for equilibration and performed 10^ measurements, 
seperated by 10 decorrelation sweeps from each other. The eigenvalues A^^^ A*^^\ A*^*^ are 

sorted in decreasing order and so the shown state A^^^ corresponds to the lightest state in the 

physical spectrum. In the r.h.s. plot, red plusses encode the contribution from to the physical 
eigenstate , green crosses stand for 0;^, blue circles for (p^x* magenta boxes for X*X- From 
the obtained results we conclude that the state A^^^ is the symmetric or anti-symmetric combination 
{(pX ± (p*X*)lV2, associated with a mass of m « 1.98, which was obtained from a fit to A^'^^ at 
An, = 1,2,3. We do not show the lightest state because it turns out to be of the form (p*(p, 
which is a combination we do not expect to be excited by the introduced iso-spin like chemical 
potential ihj, = = ji. Naively one expects that the condensation thi'eshold is related to the mass 

of the lightest state m the chemical potential couples to, by Ipc = tn, where in our case this state 
is one of the two combinations {(px <p*X*) /v^- Guided by Hand’s rule we tentatively make the 
assumption that the state with the higher multiplicity is the one with lower energy and so we expect 
the threshold of condensation to be connected to the mass of the triplet state {(px + / V^, thus 

identifying A^^^ with the symmetric combination (consistent with the r.h.s. plot in Figure 1). 

In Figure 2 we show the paiticle number densities and for the respective flavors and fhe 
corresponding susceptibilities as function of fhe chemical pofenfial. The paramefers were again sef 
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Figure 2: We show the particle densities n^ and n^j^ and the corresponding susceptibilities and as 
function of the chemical potential for j3 = 0.7, A = 1, = 5.3, = 5.7 at x N, = x 8,12. 


to j8 = 0.7, A = 1, =5.3 and = 5.7 at a spatial lattice volume of = 8^, where we used 10^ 

equilibration steps and performed 10^ measurements seperated by 10 steps. We observe a single 
threshold of condensation at /fc « 1.6 in all shown observables and also in the field expectation 
values (I^P) and {\x\^) and the plaquette expectation value {U) which are not shown here because 
of limited space. The threshold in a run at higher temperature {Nt = 8) is shifted towards smaller 
H as is expected due to the higher thermal excitation of the system. The mass of the state A^^^ and 
also the other eigenvalues of the correlation matrix C,y are plotted as vertical lines, where it can 
be seen that the critical chemical potential does not match with the mass of A^^^, as was naively 
expected. 

We will suggest an explanation for the observed threshold in the following subsection. 


5.1 Relating the condensation threshold to the lowest excitation of the system 


We define the fields Xjc and yx at position x as linear combinations of the fields (j)x and Xx 
appearing in the original action given in Section 2 


_^x + Xx _<l>x-Xx 

s/2 ' s/2 


(5.1) 


In terms of the fields x and y we can reexpress the condensate which forms at the critical chemical 
potential and which we expect, as was argued before, to be of the form (pxXx + A? 

<i>xXx + (l>xXx = M^-\yx\^ ■ ( 5 - 2 ) 


Since in the condensed phase the above expression will become large, we make the simplifying 
guess that at jic the modulus \xx\ will also grow large, while \yx\ will become very small. Then 
reexpressing the original action in terms of the fields x and y and ignoring contributions which are 
at least linear in |y| due to |y| being small, we obtain the reduced action Sc[U ,x]\ 

hm 5[[/,x,y] ^ S,[[/,x] = -j8. £ £ Re Ux.aUx+sM+^ (5-3) 
bKO x^z 

+1 [Ml k.P + A,|X.|" - £ [x: Ux,vXx+^ + X*x ■ 

^ V =1 ^ 
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Where the reduced parameters are given by 


Ml = 


Mi+Mi 




i8. = j8 


(5.4) 


Note that since Xc = ^ we expect the mass of the field x to be smaller than the masses of the original 
fields (j) and X- So in fhe condensed phase fhe sysfem is described by fhe reduced acfion Sc\U,x] 
and we expecf fhe lowesf mesonic excifafion of fhis reduced sysfem fo be relafed fo fhe fhreshold 
of condensation shown in Figure 2: Thus we expecf Ijic = mx*x, where fhe mass mx*x corresponds 
fo fhe fwo-parficle bound stafe formed by one x- and one anfi-x particle af zero chemical potenfial, 
wifh fhe dynamics of fhe field x described by fhe reduced action S'cff/jx] given in (5.3). 

We exfracf fhe mass mx*x of simulafions wifh fhe reduced acfion Sc[t/,x] af fhe reduced param- 

efers 

A/a ~\~My X 

M^= ^ =5.5 , A, = -=0.5 , j3c = j8=0.7 , (5.5) 

which are sef according fo (5.4) wifh respecf fo fhe original paramefers we used in fhe finile chem¬ 
ical pofenfial plofs in Figure 2, so we can compare fhe obfained mass fo fhe defermined fhreshold. 
The resulf for fhe mass mx*x is shown in Figure 2 as dashed vertical line labeled mx*x, where if can 
be seen fhaf fhe obfained mass does mafch fhe fhreshold very well. We lake fhis as positive lesl for 
fhe explanalion given above. 


6. Summary 

Wilhin fhe sludied model we suggesfed and verified an explanalion for fhe non-frivial relation 
befween fhe obfained He threshold of condensation and the determined /r = 0 mass spectrum. In 
further studies it would be interesting to investigate the interaction between the excited mesons in 
the condensed phase to better understand the first-order nature of the condensation transition. It 
would also be appealing to introduce a third flavor info fhe Iheory fo have a richer speefrum and as 
a consequence possibly encounfer mulliple condensafion fhresholds. 
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